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Abstract 

The anomalous action of the chiral effective theory to 0(p ) is investigated by generalizing the 
consideration in jl| with including the wave function part in the general quark propagator. It is 
found that the QCD dynamics dependence of the Wess-Zumino term is explicit and this dependence 
means that the QCD chiral symmetry must be broken dynamically in the low energy region for 
dynamically generating the Wess-Zumino term. In addition, we found that the next to leading 
order anomalous action in the momentum expansion of the chiral perturbation theory is gauge 
invariant and QCD dynamics dependent. 
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I. INTRODUCTION 



Due to the quark confinement in the low energy region, the strong interaction processes 
cannot be calculated in the frame work of conventional perturbation theory in which the 
expansion is in terms of the coupling constant g s . In this sense, to describe the low energy 
dynamics of QCD, some effective theories, based on the Weinberg's folk theorem^ that any 
quantum theory that at sufficiently low energy and large distances looks Lorentz invariant 
and satisfied the cluster decomposition principle will also at sufficiently low energy look like 
a quantum field theory, was proposed. Among all these effective theories, chiral perturba- 
tion theory^, got great triumphes in describing the hadron processes with pseudoscalar 
mesons. 

In principle, there are two wavs to obtain the chiral perturbation theory. One is the 
phenomenological construction (3, e|, in which, chiral symmetry SU(3)r x SU(3)l and its 
breaking SU(3)r x SU(3)l — > SU(3)y of QCD was considered, the pseudoscalar mesons 
emerge from the chiral symmetry breaking as Goldstone bosons and live in the coset space 
of subgroup of original chiral symmetry. The lagrangian is expressed in terms of the order 
of small quantities such as the ratio of external momenta of mesons to the chiral symmetry 
breaking scale, and in front of each independent structure term, a coefficient is given. It 
is believed that the QCD dynamics are incorporated in the unknown coefficients and these 
coefficients can be determined from the underlying QCD in principle, but in practice they 
are determined phenomenologically due to our few knowledge about the nonperturbative 
calculation. In this approach, the pseudoscalar mesons get mass through the Gell-Mann, 
Oakes and Renner relations 4] by adding the quark masses to the chiral effective theory 
which breaks the chiral symmetry explicitly. In fact, this method had been used to describe 
the hadron physics before the appearance of QCD|5[. 

Another approach is to obtain the chiral perturbation theory from QCdQ], in which, the 
relation between the phenomenological lagrangian and QCD was established. It is found 
that the coefficients of the phenomenological lagrangian can be expressed in terms of the 
Green functions of quarks so that they can be regarded as the fundamental QCD definitions 
of the coefficients of the phenomenological lagrangian. Then, unlike the phenomenological 
construction in which the coefficients were determined phenomenologically, in the QCD ap- 
proach these coefficients can be evaluated numerically from the first principle of QCD as 
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long as we have found method to compute all related Green's functions of QCD. As the 
first approximation, the explicit results were given in which can be shown further to be 
equivalent to the computation based on a gauge invariant, nonlocal, dynamical model (GND 

grangian were expressed in 

terms of the quark self-energy and the numerical results based on the quark self-energy 
determined in some models were also given, these numerical results coincide with the ex- 
perimental ones. In fact, there was another approach by using the anomaly method[9(, but 
the signs of the numerical results of coefficients L 7 and L 8 predicted in this method are 
different from the phenomenological results. In addition, it seems that all the coefficients in 
this anomaly approach is independent of the QCD dynamics. In fact, it is shown[7] that the 
anomalous contributions are cancelled by some parts of the normal contributions and the 
cancellation leaves the quark self-energy dependent contributions to the coefficients of the 
chiral lagrangian. 

What we mentioned above are mainly the normal part of the chiral lagrangian. Besides 
the normal part, there is anomalous part with odd intrinsic parity of the chiral lagrangian 
which is due to the negative parity of the pseudoscalar mesons. The anomalous part of the 
chiral lagrangian was first given in 10] by integrating the consistent condition of anomal y in 
the language of current algebra. Latter, this lagrangian was constructed geometrically llj 
with photon as an external field. Its topological meaning was studied in detail by several 
groups Like the normal part of the chiral effective theory, the quantum corrections 

do not renormalize the coefficient of the Wess-Zumino term since it is the leading order of 
the anomalous action. For the loop corrections, it is found that they are gauge invariant and 
give rise the higher order anomalyJLJ] since the Wess-Zumino term is not gauge invariant, 
the explicit form of the next to leading order were give inQ|. This does not mean that 
the nonrenormalization theorem is violated, because only the summation of all the orders 
of the effective theory is equivalent to the underlying theory. For the anomalous part of the 
chiral effective theory, as for the normal part, among all these constructions, the QCD effects 
were hidden in some constants such as the pion decay constants, so that the explicit effects 
of the QCD dynamics were also not clear. Because of all these reasons, it is interesting to 
investigate the explicit QCD dynamics dependence of the anomalous action. In our previous 
workjlj, in the fame work of GND model, we found that the anomalous terms come from 
both dynamics dependent and independent sources, after some cancellations, the Wess- 
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Zumino term is yielded and it is found that the coefficient relates to the quark self-energy. 
In addition, when the quark self-energy vanishes, the anomalous action also vanishes, this 
means that the anomalous action is a QCD dynamics dependent quantity. In fact, before 
our work, the dependence of strong dynamics of Wess-Zumino term was investigated in 
the constituent quark model by introducing a constituent quark masses Mn jl6^ and it is 
found that the Wess-Zumino term vanishes when the strong interaction was switched off by 
set Mq = 0. But in all their approximations, the hard constituent quark masses Mq will 
cause wrong bad ultraviolet behavior. The anomalous section of the effective action and its 
applications were reviewed infisl] . 

In fact, the quark propagator in GND quark model is first order approximation of the 
dynamical perturbation theory and when one switches off the external sources, only the 
quark self-energy effect is considered. Generally speaking, the inverse of quark propagator 
takes the form 

S-\k) = iA(k 2 )Jj:- B(k 2 ) (1.1) 

where B(k 2 ) arises from quark self-energy and the quark condensation which indicates the 
chiral symmetry breaking if 

{S- 1 (p), l5 } = -2 l5 B(p 2 )^0 (1.2) 

and A(k 2 ) is quark wave function part. In the high energy region, their explicit forms can 
be determined order by order explicitly through the ordinary perturbation theory, but we 
cannot get their analytic forms without approximation in the infrared region yet since we 
do not know how to deal with nonperturbation theory. 

As we mentioned before, in leading order of dynamical perturbation theory A(k 2 ) = 1, in 
addition, Schwinger- Dyson equation tells us that A{k 2 ) = 1 only happens in Landau gauge 
with bare gluon-quark vertex and bare gluon propagator in the kernel of the integration 
equation. So that, to respect real original QCD effects, we have to include the A(k 2 ) effect. 
In the previous work Q and jjj, it is for simplicity of computation, we take the approxima- 
tion A(k 2 ) = 1. The main purpose of present paper is to overcome this shortcoming and 
investigate the effect of A(k 2 ). 

In literatures, Z(k 2 ) = 1/A(k 2 ) and Effc 2 ) = B[k 2 ) / A{k 2 ) are called the quark dressing 
function and mass function respectively[23[]. In the language of renormalization, the quark 



propagator (jl.lj) depends on the renormalization point /i, i.e. 

S-\n, k) = iA(n, k 2 )^ - B{n, k 2 ) 
and in the standard momentum subtraction scheme 



(1.3) 



Z(fi,fi 2 ) = 1 

E(fi 2 ) = B(fi,fi 2 )/A^,fi 2 ) 



m(/x) 



(1.4) 



(1.5) 



And the scale [i characterized the mode of chiral symmetry breaking. As was shown in|l9j|. 
when the scale is large enough which makes QCD in the perturbation region, quantum cor- 
rections cannot break the chiral symmetry dynamically if the current quark mass which 
breaks the chiral symmetry explicitly was not introduced. While in the low energy, the 
chiral symmetry is broken due to the quark condensation and the mixed quark-gluon con- 
densation. In case that the current quark mass is introduced, the chiral symmetry will be 
broken explicitly due to the small light quark masses when scale /i is large, but there is 
no Wess-Zumino anomaly generated. When [i becomes so small that the theory is in the 
nonperturbative region, the dynamical chiral symmetry breaking term becomes larger than 
current quark mass and the the chiral symmetry is broken dynamically. Considering this, 
at high energy region, the renormalized quark propagator can be related to the bare one 
through the wave function renormalization constant Z 2 through 



where c is a small constant which is induced in a regularization such as e in dimensional 
regularization and describes the pole behavior of Z%. Since the physical or the renormalized 
propagator S^/i, k) is finite, the bare propagator should be negative power of c which is 
needed to cancel the infinities arise from the loop integral. Then, for a sufficiently small c, 
the relations between the wave function renormalization Z 2 and the quark dressing function 
Z at different renormalization point can be yielded as 



which means that mass function must be renormalization point independent and a change of 
the renormalization point is just an overall rescaling of Z(fi, k 2 ) by a momentum independent 





Z 2 (fi2,c) Z(^c) 



(1.7) 




(1.8) 
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constant. Then, the measure of the nonperturbative physics is the deviation of Z(fi, k 2 ) 
from 1 and the difference of S(/c 2 ) from the renormalizaed quark mass m(p) j^. In this 
paper we will not distinguish the concepts of quark dressing function and the wave function 
renormalization constant and call A(k 2 ) as the wave function renormalization constant. 

Formally, the functions A(k 2 ) and B(k 2 ) are constrained by the Schwinger- Dyson equa- 
tion. Although Schwinger- Dyson equation may give a strong constraint on the dress function 
A(k 2 ) and quark self-energy B(k 2 ), it cannot be solved exactly since it is a recurrent inte- 
gral equation. Based on some approximations, special forms of A(k 2 ) and B(k 2 ) were got 



from the Schwinger-Dyson equation in Landau gaugejTJ, 20], axial gauge |2!| and covariant 



gauge 



221. 



In this paper, we will study the anomalous action based on the general quark propagator 
(II. lj) without discussing the implications of the dynamical chiral symmetry breaking which 
have been discussed in|25( and also the amplitude of the scale /i, we only consider the theory 
in the nonperturbative region. For the convenience of our discussion, we will introduce 
the vector, axial-vector, scalar and pseudoscalar sources like what was done in the normal 
section of the chiral effective lagrangian! 3]. The anomalous action we will investigate is the 

nn 

part of the action proportional to odd number of Levi-Civita tensor e^ aj g[18, 26]. 

This paper was organized as follows: In section |Hj we will discuss the meaning of chiral 
symmetry by considering the relation between the chiral effective lagrangian and the effective 
action we are interested in and give chiral invariant general quark propagator. In section HTT| 
we shall calculate the Wess-Zumino term from the general quark propagator in section [H] 
and investigate the effects of QCD dynamics in the anomalous section of the chiral effective 
lagrangian. Our conclusions and outlooks are presented in Section ITVl 



II. VARIATIONS OF THE EXTERNAL SOURCES UNDER CHIRAL TRANS- 
FORMATION. 

In this section, we shall analyze some transformation properties of the external sources 
under chiral rotation that we will use in the following calculations. The explicit action which 
is proportional to the inverse of full propagator (jl.ljl is presented when we switch off the 
external sources by comparing the Gasser-Leutwyler lagrangian [^. 

It is found that the introduction of some external sources in the study of the chiral 
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perturbation is convenient. Following ref. in four dimensional Euclidean space-time, we 
consider pure quark kinetic part of bare QCD action in the presence of external scalar, 
pseudoscalar, vector and axial-vector sources 



S[ip,ip,J] = J d 4 xipDiP 

= J d 4 x[ip^ + ^J(x)tlj] (2.1) 

where the external sources are defined as 

J(x) = —if{x) — id(x)j 5 — s(x) + ip(x)^ (2.2) 

and quark mass matrix can be extracted from the scalar external source s. In chiral pertur- 
bation theory, the power counting of the vector t> M and axial- vector a M sources are both 0(p) 
while scalar s and pseudoscalar p sources are both 0(p 2 ). For the following convenience, we 
rewrite the Dirac operator D as 

D = f-s + ip^ 5 , f=i t *V ll V M = dp - iv„ - 10^75 = -V+ (2.3) 

Formally, the generating functional of QCD is 

Z[J\ = J Vil)V$V^V^VA^^ x[c ^^^> A ^ J ^ (2.4) 

where C(ip, ijj, *ff, *ff, A^) is lagrangian of QCD with ip, \l/ and are light, heavy quarks and 
gluon fields respectively. By integrating out the quark and gluon fields and integrating in 
the pseudoscalar mesons, one can formally rewrite the QCD generating functional as 

Z[J\ = [ VUe s ^ [u > J] (2.5) 



where Sql[U, J] is the chiral effective lagrangian of pseudoscalar mesons with external 
sources. In principle, Sql[U, J] consists normal part and anomalous part due to the in- 
trinsic parity of pseudoscalar mesons, that is 

£>Gl{U, J] = S norma \[U, J] + SznomalyiU, J] (2.6) 



the normal part S norma i[U, J] w as g iven in to 0(p 4 ) and the leading anomalous part 
^anomaly [Uj J ] can be found infl [J and next to leading part was given inQ. 



Then, the action (|2.1|) which we are interested relates the generating functional through 

Z[J\ = f VUVi)Vi)e Scs[ ^^ J] (2.7) 



By comparing ()2.5|) and (|2.7jh we yield the relation between the chiral effective action 
()2.5j) and the effective action ()2.7j) 



Because the massless QCD is invariant under chiral transformation, the action 
S e s[ip,ip, U, J] is required to be invariant under the following local SU(3) R x SU(3)l chiral 
transformation 

$(x) -> ij'(x) = [V L (x)P L + V R (x)P R ]i/;(x) (2.9) 
J{x) -> J'(x) = [Vl(x)Pr + V^(x)P L ][J(x) + ^][Vj(x)Pi + vJ(x)Pr] (2.10) 
U'(x) = V R (x)U(x)Vl(x) (2.11) 

where Vl{x) and Vr(o;) are the left and right chiral rotation matrices and P R) l are the 
project operators \{1 ± 75). 

Mathematically, the pseudoscalar field U(x) has the decomposition U(x) = Q 2 (x), under 
local chiral transformation SU(3) R x SU(3)l, Q(x) transforms as 

Q(x) -> fi'(x) = h\x)Q{x)Vl{x) = V R {x)Q{x)h{x) (2.12) 

and the gauge group h(x) can be determined by Vl, V r and Q(x). The symmetry /i(x) is the 
so called hidden local symmetry, its gauge bosons can be identified with the vector mesons 
through the so called hidden local symmetry method (22], 13- 

In terms of the nonlinear field Q(x), we can introduce the rotated external fields 

= -ii'n(x) - - s a (x) + ipn{x)j 5 (2.13) 

Replacing all the external fields with the rotated ones and defining 

D a = %-8n + i Pa jB, V^ = d^ -iv^-ia^ 5 (2.14) 

then, under chiral transformation SU(3) R x SU(3)l, one can easily prove the transformations 

Jq(x) -> J' a (x) = h\x)[J n (x) +$ x )h(x) 

D Q -> D' Q = h\x)Do,h{x) (2.15) 
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From above transformations, one can conclude that the local chiral transformation can be 
realized through hidden local symmetry, once the theory keeps the hidden local symmetry, 
it preserves the chiral symmetry automatically. That is, by introducing the Q field, the 
chiral symmetry can be reflected through the hidden local symmetry In fact, even in the 
presence of massive quark, by replacing the external fields with the rotated ones, the local 
chiral symmetry can be preserved 3^ |. 

For the convenience in our discussion of the Wess-Zumino term, we extend Q(x) into five 
dimensional space [33J] with 



Q(t, x) = e itxa * a{x) Q(x) = x) £1(0, x) 



(2.16) 



and define the t dependent rotated sources 

Jn(t,x) = -if n(t,x) - i$ a{t,x)j 5 - s Q (t,x) + ipn{t,x)^ 5 

= [n(t,x)P R + tf(t,x)P L ] [J(x)+p x ] [Q(t,x)P R + Q\t,x)P L ] (2.17) 



we can verify 



sn(t,x) 



-[Q(t,x)[s(x) -ip(x)]Q(t,x) + &(t,x)[s(x) + ip(x)}Vt\t,x)} (2.18) 

~[tl(t,x)[s(x) -ip(x)]tt(t,x) - &(t,x)[s(x) + ip(x)p\t,x)} (2.19) 

i[fi f (t,x)[f M (x) +a M (x) + id^}Q{t,x) 

+^(t,x)[v IM (x)-a fl (x)+id fl ]^(t,x)} (2.20) 
^[Q j (t,x)[v^(x) + a^(x) + idjft(t,x) 

-n(t,x)[ VlM (x) -a p {x) + id^\t,x)} (2.21) 



Now consider infinitesimal transformation at parameter space t — > t + St, correspondingly 

Q(t, x) -> Q(t + St, x) = (l + iSt\ a n a (x))tt(t, x) = n(t, x)(l + iSt\ a n a (x)) (2.22) 
Q\t,x) -> rt(t + 6t,x) = (l-i6t\ a n a (x))n\t,x) = rt(t,x)(l-i6t\ a Tr tt (x)) (2.23) 

we can easily prove 



SJn(t,x) 
8a^(t,x) 



{iSt\ a Ti a (x)j 5 ,Jn(t,x) + $ x } 



St 

"2 



su 
It 



U\v^{t,x)] 



iSt SU 



tf(t,x),a&(t,x)} 



(2.24) 
(2.25) 

(2.26) 



And in addition 

5%(t,x) = -|{[AV(x),a^(t,x)]V^(t,x) + V^(t,x)[A a 7r a (a;),a^(t,a;)]} (2.27) 
For the dimensional extended Dirac operator Dn(t,x) 

D n (t, x) = - suit, x) + ip u (t, x) 75 = p + J Q (t, x) (2.28) 

we have 

6D a {t,x) = {iSt\ a 7r a l5 ,D n (t,x)} = \{ W ^l ^ U\t, x) l5 5t, D a {t, x)} (2.29) 

2 dt 



Equivalent ly 

dD n (t,x) 



dt 2 L dt 



= -\{U\t,x)^fl l5 ,DUt,x)} 



(2.31) 



By now, the transformation properties of the external sources were derived from the 
classical lagrangian (|2.1j) . we want to include the QCD dynamics effects next. These ef- 
fects can be calculated perturbatively at high energy region, however, at low energy region, 
the quark-gluon coupling g s is strong and it is not suitable to regard it as an expansion 
parameter. In this sense, we do not know how to get the forms of the wave function renor- 
malization constant A(k 2 ) and the quark self-energy B{k 2 ) although formally they satisfy 
the Schwinger-Dyson equation. Since we have the conclusion that the Wess-Zumino term 
is an intrinsic property of the non-Abelian gauge theory both phenomenologicallyjnj and 
geometrically [ill , we naively want to know whether we can get the Wess-Zumino term when 
the QCD non-Abelian gauge interaction was included in the action (jl.lj) and explore the 
contributions of the wave function renormalization constant A(k 2 ) and self-energy B(k 2 ). 

To exhibits the form of the action include QCD non-Abelian gauge interaction and its 
dynamical effects, following the method given in ref. [8J], let us turn to the relation ()2.8|) 
between the Gasser-Leutwyler action S'gl and the effective action S e ff- By using the rotated 
basis, this relation can be rewritten as 

,S eS [^,U,J] 
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where S[ip, 4>, J] was defined by (|2.1|) and Af = det[(j) + J}. The anomaly due to the variation 
of the integration measure was cancelled between the denominator and the numerator. The 
relation (J2.32)) means that, formally, we can write S e g as 

Sfsi^a^aAtM = S[ipn,iJQ,Jn] + S- mt [ip n , J n ] (2.33) 

where S[ipn, ipri, Jn] is the action defined in (j2.1|) by substituting relative quantities with 
the rotated ones. S^t are the terms represent effects of color gauge interaction which are 
proportional to a s at the leading order, it includes the fermion self-energy term and the wave 
function renormalization constant which are caused by integrating out the heavy fermion 
and gluon degree of freedom in the underling QCD and integrating in the local pseudoscalar 
fields U(x). We formally write the interaction part Si nt as 

SUVta, M ~ / d 4 xMx){[A(d 2 ) - + Jn) + B(d 2 )}M%) (2.34) 

Where we have taken the minimal extension by including in effect of quark wave function 
renormalization corrections into the Si nt and 

J n = -if n(t,x) -if, a{t,x)j 5 (2.35) 

It should be noticed that (j2.34j) is the most general form when the dynamical effects are 
independent of external sources. In the presence of external sources, the general effective 
action (|2.34|) due to the color gauge interaction can be decomposed in terms of its general 
spinor structure which will introduce many more other functions besides the present A(d 2 ) 
and B(d 2 ), these extra functions will cause much more computation difficulties and we will 
leave the discussion of them in future. Now as first step beyond original quark self-energy 
and reduce the difficulties of the computation, we only consider one more extra function: 
quark wave function renormalization function A(d 2 ) but leaves the discussion of most general 
propagator elsewhere. In (J2.34j) . we select the dynamical dependence coefficients of the 
differential operator and the vector source are the same because we need the variety of 
vector sources to compensate the variety induced by the differential operator under the 
local chiral transformation, the coefficients of the vector source and the axial-vector source 
are also the same as we considered the combinations 

Vg„ = \K + a%] (2.36) 

Vt = \ltf-4\ ( 2 - 37 ) 
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and the Parity invariance of QCD induced property 

P ■ V£, - V£, (2.38) 

Combining this with (|2.33|) . we yield the effective action as 

S eS [ipa,ipa, 1, J n ] ~ / d 4 £ ^(x){A(9 2 )[^ + J n ] + s n - ip 7fi + B(d 2 )}M^) (2-39) 

Substituting ()2.39|) into ()2.32j) . we get the full expression of Sql as 

S GL [U,J] ~ lndet{A(9 2 )[^ + J n ] + ^-^75 + 5(9 2 )} 

- In det + J n ] + In det + J] (2.40) 

It should be noticed that there is no factors A(d 2 ) and B(d 2 ) in the last two determinants, 
this is because these two factors originate from the color interaction and the last two deter- 
minants are from the free quark lagrangian. 

As we discussed before, the local chiral symmetry can be reflected by the hidden local 
symmetry through the introduction of the rotated sources, but it is obvious that the action 
(j2.40j) cannot preserve the chiral symmetry because 

h\x)A{-d 2 )h{x) = A[-h){x)d 2 h{x)} = A[-[d^ + h\x)d lx h{x)f] (2.41) 

Similarly for B(—d 2 ). This means that the chiral transformation always induces an extra 
term in A{— d 2 ) which makes the theory change under chiral transformation. To compensate 
this extra term we consider the operator A(— Vq) in stead of A(— d 2 ) in the action with 
Vq = <9 M — Wq. Because, under local chiral transformation 

V^^V^ = h\xMh(x) (2.42) 

that is, the variation of the vector source compensates the variation induce by the action of 
the differential operator on the local phase factor, we have 

A(-V 2 n ) - A(-V 2 ) = h\x)A(-n)Hx) (2.43) 

which indicates that A(— Vq) transforms homogenously under chiral transformation. Using 
the same analysis on B(—d 2 ), we finally get the local chiral invariant action as 

Sgl[U,J] ~ lndet{A(-V 2 n )i$ + j n ] + sn-^Pnl 5 + B(-V 2 n )} 

- In det $ + J a ] + In det [0 + J] (2.44) 
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Now let's make some comments on this action: When we do not consider the wave function 
renormalization, i.e., set A = 1, we arrive at the GND quark model and its anomaly had been 
analyzed in[lj. If the QCD dynamical effects were switched off, that is, A = 1 and B = 0, 
the first two terms of this action cancelled each other and the left term is independent of 
the pseudoscalar fields, then one cannot yield the Wess-Zumino term in this case. As shown 
inQ, the scalar source s and pseudoscalar source p do not contribute to the leading order 
anomaly, so that we will not consider these two sources in the following, therefore, our final 
action can be written as 

S GL [U, J] = In det [D' n ] - In det + J n ] + In det [$ + J] (2.45) 

with 

D'n = A(-V 2 n )% - B(-V 2 n ) (2.46) 



III. CALCULATION OF ANOMALOUS EFFECTIVE LAGRANGIAN WITH 
THE GENERALIZED FERMION PROPAGATOR 

In this section, based on the action (|2.45|) . we will calculate the anomalous section of 
chiral effective action to the leading order, 0(p A ). Since the difference of (I2.45|) with that 
in our previous work Q| is quark wave function renormalization function A, the questions 
are interested become that: is there any changes due to arbitrary A for the leading order 
anomalous action? Does this anomalous action has its QCD dynamics dependence? Is there 
higher order anomalous action and where is it from? What is the meaning of higher order 
anomalous action? 

Our consideration is: If there is anomalous action, it should be evaluated through the 
direct calculation of the action (|2.45j) . To investigate the QCD dynamics effects in the 
anomalous action, one can set A(k 2 ) = 1 and B(k 2 ) = in the action ()2.45|) which mean 
that the strong interaction effects were switched off. In addition, if there is anomalous action 
by the total calculation of 1)2.45)1 . this anomalous action can be evaluated through the terms 
in the r.h.s. of 1)2.45)1 . so we need to calculate the anomalous action from every term in the 
r.h.s.of 1)2.45)1 . Since the first term in the r.h.s. of (j2.45j) is QCD dynamics dependent, we 
expect the exact form of the Wess-Zumino term and the Bardeen anomaly will impose some 
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constraints on the behavior of the strong dynamics dependent function A(k 2 ) and B(k 2 ) 
even though this constraint may be rude. 



A. Direct Calculation of Anomalous Action Based on the Effective Action (|2.45j) . 



To investigate the anomalous action in the effective action ()2.45|) . we make an inverse 
rotation of (|2.13|) which makes the last two terms of ()2.45|) cancel each other and the pseu- 
doscalar fields in the second term rotate to the first term, that is 



where 



S GL [U,J] = lndet{A(-V 2 n W-B(-Vl)} 



A(-Vl) = p(t,x)P R + Q(t,x)P L }A(-V 2 n )[Q(t,x)P R + Q\t,x)P L ] 



B(-wl) = p(t,x)p R + n(t,x)p L ]B(-vl)p(t,x)p R + n(t,x)p L ] 

Considering identities IndetA = TrlnA and 5TrlnA = Tr5AA~ l , we have 



S[U, J] - S[l, J) 



dm 



dt nK ' 



where the trace Tr is around the color, configure, flavor and Lorentz space and 

D' n (t) = A(-Vl)f-B(-V 2 n ) 



(3.1) 

(3.2) 
(3.3) 

(3.4) 

(3.5) 



It should be noticed that we introduced a term S[l, J] in ()3.4|) which would not change our 
conclusion on the Wess-Zumino term since it is independent of the pseudoscalar fields. 
Since the anomalous effective lagrangian is proportional to the odd number of Levi-Civita 



tensor e^ vaP 



261 ] . i.e., odd number of 75 in the trace of Lorents matrices, in the following we 



only concentrate on this kind of terms. After the trace around configure space, we have the 
anomalous action as 



S[U,J]-S[1,J] 



dt id x 



d 4 k 



rtr. 



c,f,l 



3D', 



ot 



(3.6) 



/o J J (2tt) 4 

the subindex e means we are only interested in the terms consist odd number of Levi-Civita 
tensors and d — > ik + d indicates that after the trace on the configure space, the differential 
operator d should be replaced by ik + d, that is 



D'na-Mk+a = A[-(y n + ik) 2 }\y+ ik 1 } - S n [-(Vn + ik? 



(3.7) 
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The next problem we have to solve is how to deal with the function A[— (Vq + ik) 2 } and 
B[— (Vq + ik) 2 ]. To illustrate this problem, let us focus on A[— (Vn + ik) 2 ], to deal with this 
function, one has to deal with function A[— (Vq + ik) 2 ] at first. To make its power counting 
explicit, we expand it around k 2 , then using Taylor expansion, we have 

A[-(tk + Vn) 2 ] = A[k 2 -2tk-V n -V n ] 

oo / i \n 

= A(k 2 ) + V A (n \k 2 )^-L(2ik ■ V n + V 2 n ) n 

»=i n! 
= A(k 2 ) + A(k„V^) (3.8) 

where A is a function of external source and its chiral counting is equal or higher than 0{p), 
explicitly, it is 

Vg) = £ A^(k 2 ) [ —L(2tk ■ V Q + Vl) n (3.9) 

n=l n - 

So that 

A[-(Vn + ik) 2 } = A(k 2 ) + A(k„ V£) (3.10) 

with 

A(k^,V& = p(t,x)P R + Q(t,x)P L ]A(k II ,V^[Q(t,x)P R + Q\t,x)P L ] (3.11) 
Similarly, we can write 

B[-(ik + V n ) 2 ] = B(k 2 ) + B(k„ V£) (3.12) 



with 



Then 



with 



Bfa, Vg) = £ 5< B >(fc 2 )^-/-(2ifc • + V 2 )" (3.13) 

71=1 n - 



B[-(V n + ^) 2 ] = B(k 2 ) + V£) (3.14) 



B(k 2 ) = B(k 2 )[U\t,x)P R + U{t,x)P L ] (3.15) 
B(k^V& = p(t,x)P R + Q(t,x)P L ]B(k,,V^)p(t,x)P R + Q(t,x)P L ] (3.16) 
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Then, we can formally write 



D'nit) = A(k 2 )tlf - B(k 2 )[U^(t,x)P R + U(t,x)P L ] + £[V£,. 



(3.17) 



with 



£[V£,a£] = A(k^)0 + {A(k 2 ) + A(k^)}f -B(k^M 



(3.18) 



So that 



A{k 2 )i^ - B(k 2 )[U\t,x)P R + U(t,x)P L ]+S[V^,a, 
{A(k 2 )t^ - B(k 2 )[U\t,x)P R + U(t,x)P L ]}- 1 
x{l+^[A; /1 ,V^[A(A; 2 )^- J B(A; 2 )[f/t( t)X )p i? + f/ ( t)X )p L ]]-i}-i 
-A(k 2 )i1jk - B(k 2 )[U(t,x)P R + W(t,x)P L ] 



x 



£(-i) m {£[v£ 



A 2 (k 2 )k 2 + B 2 (k 2 ) 

A(k 2 )i^ - B{k 2 )[U{t, x)P R + W(t, x)P L ] 



AM 



^5 "n 



} m (3.19) 



A 2 (k 2 )k 2 + B 2 (k 2 ) 

By now, the calculation becomes straight forward although complex. Explicit calculation 
gives 



S[U,J]-S[1,J] 



N C C 

■ 4g7r2 W 



dt I d 4 x 



xtr, 



U\t, x) dU ^^ L^t, x)L u (t, x)L a (t, x)Lp(t, x) 



(3.20) 



where N c is the number of color, L^{t,x) = W(t,x)d^U(t,x) and 



C = — I d k 



4A 4 (k 2 )B%k 2 ) 



[A 2 (k 2 )k 2 + B 2 (k 2 )] 5 [A 2 (k 2 )k 2 + B 2 (k 2 )} 5 
A 2 {k 2 )B\k 2 )k 2 d B 2 {k 2 ) 



8A 4 (k 2 )B 5 (k 2 )B'(k 2 )k 2 8A 3 (k 2 )A'(k 2 )B 6 (k 2 )k 2 



— 12 / dk 2 

Jo [A 2 {k 2 )k 2 + B 2 {k 2 )} 3 dk 2 [A 2 {k 2 )k 2 + B 2 {k 2 )} 



[A 2 (k 2 )k 2 + B 2 (k 2 )f 

(3.21) 



This constant shows the QCD dynamics dependence of the Wess-Zumino term 



Comparing this conclusion with that given in[l| , it is seen that when we set A(q 2 ) = 1 
which means the leading order of dynamical perturbation theory, the present result is the 
same as that ofjlj]. From the present expression (|3.2(J|) . it is clear that the Wess-Zumino term 
arrives from the action (J2.45|) which depends on the dynamics of QCD, when the dynamics 
of QCD is switched off by taking A(k 2 ) = 1 and B(k 2 ) = 0, the coefficient C also vanishes, 
this means that the Wess-Zumino term is an intrinsic property of QCD and is a result of 
the strong interaction. Note that the QCD dynamics dependence of Wess-Zumino term has 
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also been investigated in 29| based on the GCM model, but our result here has a factor 



y. 



k 2 A(k 2 )/[k 2 A(k 2 ) + B 2 (k 2 )} difference with the one in 

In the case of dynamical chiral symmetry breaking, i.e. B(k 2 ) ^ 0, the pseudoscalar 
mesons emerge as the Goldstone bosons and there will be anomalous processes which can 
be described by the Wess-Zumino term^oj]. In this case, one should have C = 1 to recover 
the exact coefficient of Wess-Zumino term. We find that as long as the functions A(k 2 ) and 
B{k 2 ) behave as 



t 



A 2 (k 2 )k 2 
B 2 (k 2 ) 



oo 




k 2 
k 2 



oo 




(3.22) 



The coefficient C can explicitly be integrated out 

t d 



C = -12 



dt- 



1 



1 



(3.23) 



io [1 + t] 3 dt [1 + t] 

which does coincides with the exact coefficient of Wess-Zumino term. 

Although relation ()3.22)1 seems like a constraint on the behavior of the function 
B(k 2 )/A(k 2 ), while it is very crude. From the explicit calculations given in the literatures, 
we can extract the ratio 



B{k 2 )/A{k 2 ) 
k 2 

B(k 2 )/A(k 2 ) 
k 2 

B(k 2 )/A(k 2 ) 
k 2 



4m|, -P. 



12 



k 4 



3 1 



7 



33 -2n 



t\-± + i^t+0(t*)} ' * u > 
A- 2 



[l-iAt-i+OCt- 2 )] 1 / 2 
rTTT] 



OO. 



ci$(2,f 



-fc 2 

2' f3M 2 



+ C 2 ( 



(3M 2 



-3/2$Cl ^1. 



/3M 2 



fi^M^ ! (^)- 3/2 'f(i- 



fc 2 



-; :\ i - 



0; 



oo. 



(3.24) 
(3.25) 

(3.26) 



where ()3.24|) was given in|l7| by using the Landau gauge, eq. (|3.25jl is from the conclu- 
sion given i„Q ba S ed on the covariant gauge and eq-ES was the result from the axial- 



gauge 



2l\ . In (|3.24|1 . n and are the number of flavors and constant dynamical quark 



mass respectively. In (13.37)1 . A is a constant due to the ghost self-energy at zero point and 
t oc k 2 . $ is the confluent hypergeometric function, (3 oc a s and M is the renormalization 
point in ()3.26|) . After some basic algebra, it is seen that the above three results (|3.24H3.2'B|) 
all satisfy the condition eq. ()3.22|) . 

To rewrite the four dimensional Wess-Zumino term (j3.20|) in five dimensional space-time 
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HI 



with four dimension space-time boundary, we use the same trick as 
§£tr f {Li(t, x)Lj(t, x)L k (t, x)L x {t, x)L m (t, x)}e ijklm 

= 5^tr f {W(t, x)^|^L,(t, x)Lfa x)L k (t, x)Li(t, x)}e ijklm (3.27) 
where is a totally antisymmetric tensor. This relation yields 

J tfE <J ' Wm tr/{Li(l, x)Lj(l, x)L k (l, x)U{l, x)L m (l, x)} 

= 5^ d 4 x J 1 dte^ af3 ti f {U^t, x) ^ L„{t, x)L v {t, x)L a (t, x)L p (t, x)\3.28) 

From this and using C = 1 we can get the standard Wess-Zumino term[lo| 

N r 

r ~[^] = / dE^ m tr / {L 4 (l,x)L J (l,x)L fc (l,x)L / (l,x)L m (l,x)} (3.29) 

Z4U7T Z JQ 

In conclusion, it is seen that the Wess-Zumino term is an intrinsic property of QCD and 
depends on the QCD dynamics. All the QCD dynamics effects in the Wess-Zumino term 
can be collected into the coefficient C which was given by 1)3.21)1 and, when the coefficient 
approaches to unit, one arrives at the exact Wess-Zumino term. When one switches off the 
QCD dynamics, the coefficient will vanish, which means a vanishing anomalous action and 
no anomalous pseudoscalar processes occur. 



B. Anomalous Action from Each Term of the r.h.s. of Effective Action (12.45(1 . 



In above subsection, we have proved that there is leading order anomaly in the action 
(12.45(1 . now let's investigate the imaginary part of the action originates from every term of 

(JUS). 

Let us look at the anomalous action from from the first term of (|2.45j) . By using the 
same method as that used in above subsection, explicit calculation shows that the anomalous 
action T^[ D due to the QCD dynamics of the first term of ()2.45j) is 



Y l D = MC^ap J dt J d 4 x^tr f 



+4iC 2 e flual3 J dt J d^x^tTf 
+4iCg,ep Va p j dt J d 4 xj^tr/ 



(3.30) 
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where the index Q means that all the external sources are rotated sources and depend on 
the chiral angle Q. And the definitions of the constants are defined as 

d 4 k A 2 (k 2 )B 4 (k 2 ) d B 2 (k 2 ) 

(3.31) 
(3.32) 
(3.33) 



( | ' {2n) 4 [A 2 (k 2 )k 2 + B 2 {k 2 )fdk 2 [A 2 (k 2 )k 2 + B 2 (k 2 )] 
r d 4 k 2A 4 (k 2 )B 2 (k 2 )k 2 d B 2 {k 2 ) 

2 ~ I (2vr) 4 [A 2 {k 2 )k 2 + B 2 (k 2 )f dk 2 [A 2 {k 2 )k 2 + B 2 {k 2 )} 

r d 4 k A 6 (k 2 )k 4 d B(k 2 ) 

3 ~ J (2vr) 4 [A 2 (k 2 )k 2 + B 2 (k 2 )f~d^[A 2 (k 2 )k 2 + B 2 (k 2 )] 

Using (|2.25jl and (|2.2fijl . we can prove relations 



d 







Vfj,V u V a ap 



I QJJ 

2 e nv a ptTf{—U '[2a p V v V a ap + V,jV v a a ap - a^V u a a V 



-V^Vaap + a^aJVaVp + 2V ^V U V a V ^n} 

QJJ 

■e^ va pti f{—U ] [2a^a u a a af } + 2V ^a u a a V 'p + c 



(3.34) 
(3.35) 



we get the anomalous action with form 



T 1D = 2N c C 1 e /J , ua/ 3 dt J d 4 xtr f 



dU 
~dt 



x{2V M V J/ V Q V / 3 + 2a tl a v V a Vp - 2V fl a u V a ap + 2V ^a u a a V p 
+2a ft V v V a ap - 2a ll V v a a Vp + 2V^V u a a ap + 2a fJ ,a u a a ap}u 

pi r vQJJ 

+2N c C 2 e imxP J dt J d 4 xtr f — IP 

x{4V^V I ,V a V / 3 + 2a ll a v V a Vp ~ 2V fl a u V a a p + 4a M V i ,V a a /3 
-2a.nVva.aVi3 + 2V fJ ,V u a a a p } n 



+2N c C 3 e^ ua/3 dt J d 4 xtr f 
x{2V^V u V a V/3 - 2V^a u a a Vp + 2a^V u V a a (3 - 2a^a v a a a f3 } n 



By using the behavior ([3.22)1 . we can exactly get 

Ci = Co = — C"\ = t: 

3 12 16tt 2 

With this equation, Anomalous action with external gauge fields can be yielded 
1 n 



(3.36) 



(3.37) 



D 



32vr 2 



J dttTf 



X 



[2% 4 8i 4 



(3.38) 



19 



where 



Vji = drf - d v v* - i[v% «?] (3.39) 
d»c% = d^-i[v%a^\ (3.40) 

which has the same form as the gauge anomaly by replacing U with the gauge transformation 

From the above relations such as ()3.30|) and ()3.36)1 . we see that the quark dynamics plays 
an important role in the form of anomalous action, if there is no quark condensation, there 
is no anomalous action as mentioned injlj. The role of quark wave function renormalization 
function A(k 2 ) is passive which play the key role only with quark self energy B(k 2 ) in the 
combination of S(/c 2 ) = B(k 2 )/A(k 2 ). In addition, the conclusion ()3.38|) means that the first 
part of ()2.45|) contributes to the imaginary part of the action ([2.45)1 . 

We now investigate how to get the Wess-Zumino termQ] from the the anomalous action 
(13.36)1 . The Wess-Zumino term is the pure pseudoscalar mesons lagrangian, so that we 
should switch off the external fields by setting J = 0. In this case, there are integrable 
conditions 

YST = iK,a v n] (3-41) 
d»a" a = <fa& (3.42) 
a£ = -n\t,x)[d^U(t,x)]^(t,x) (3.43) 

Then, substitute the above three relations to the non-Abelian anomaly (|3.38jl . we get 



fiua/3 

xti f {W{t, x) dU ^ X > L^{t, x)L u {t, x)L a {t, x)Lp{t, x)} (3.44) 

where L^(t,x) = W(t,x)d tl U(t,x). 

By now, can we say that the leading order anomalous action of ()2.45j) comes from its first 



part? As argued in 10], the anomalous of term of the effective action is not chiral invariant 
which seems conflict with the argument that the first term of (|2.45J1 is chiral invariant, 
does this means that our calculation is wrong? The answer is NO! This is because, we 
only considered the imaginary part of the first term of ()2.45)1 that depends on the QCD 
dynamical effects, if one switch off these effects in the first term of (|2.45jl through A(k 2 ) = 1 
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and B(k 2 ) = 0, one can easily find that there is another term which dependent on the 
rotated sources which will give rise the anomaly obviously but with a negative sign. Its 
contribution r^ can be calculated by using the Fujikawa's path integral method which is 
shown below. Then, one totally has 

rr = rr i0 + rr )I? (3.45) 

So that, totally, the leading order imaginary section of the action induced by the QCD 
dynamical effects in the first term of (|2.45|) cancels that induced by the rotated sources and 
leaves the next to leading order anomalous action, and, the next to leading order anomalous 
action is chiral invariant. 

Next, let us investigate the anomalous T^" induce d by the last two terms of f)2.45|) . This 
anomaly can be evaluated by the Fujikawa's method 34|.,i.e., 

T" = — Tr \n[p+ J n ] + Tr \n[p + J] 

= -/W{^p + j n (t)]} 

= - lim /Wr{^f/t 75exp[ I£±M ]} (3 . 46) 

A^oo Jo Ot A z 

After the trace on configure space, expanding the exponent, taking the limit A — > oo, 
keeping the 0(p 4 ) terms and the terms proportional to e^ ua p, we finally get the result 



T 2 = -Tr h# + J n ] + Tr h# + J] 

N c _ r 1 ,. r ,„ (dU. 

2 € vvaf3 



327T 2 

X 



4 2i 8i 4 1 ~| 

L 3 . — .3 3 3 JnJ 

which is the famous Bardeen anomaly j3J|. After using the integrability condition (j3.41H3.43*|) . 

we can get the Wess-Zumino term as mentioned above. 

In summary, from the r.h.s. of ()2.45|) . the anomalous action is 

r- = rv, + + r 2 - (3.48) 

with r^ D = r 2 = -rr )0 . 

IV. CONCLUSIONS AND OUTLOOKS 

In conclusion, start from the general quark propagator, after chiral rotation, we get the 
anomalous section of the chiral effective action with external sources. For which we can 
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yield the Wess-Zumino term after switch off the external fields. 

The anomalous action we obtained depends on the QCD dynamics closely, when we switch 
off the QCD dynamics, the Wess-Zumino term will vanish. That means that the anomalous 
processes described by the Wess-Zumino term is a QCD processes. To investigate which 
term the anomaly arise from, we explicitly calculated the anomaly induced by each term of 
action ([2.45)1 . We found that the first term of ()2.45|) can induce QCD dynamics dependent 
and independent anomalous action but with a sign difference then the leading order anomaly 
vanishes, this means that the first term of ()2.45|) only contribute to the next to leading order 
anomaly and the next to leading anomaly is chiral invariant. We also found that the last 
two terms of (|2.45j) also contribute to anomaly and with the same sign of Wess-Zumino term 
when we switch off the external sources. 

Then, similar to the conclusion given injlj], one observes two kinds of alternative cancel- 
lations: One cancellation is regarded as that the anomalous action arising from the QCD 
dynamical independent section of the first term of ()2.45j) is cancelled by the anomaly due 
to the dynamical dependent section of the first term of ()2.45j) and leaves the anomalous 
action from the second term of ()2.45|) . In this cancellation, the anomalous action is found to 
be QCD dynamics independent. An alternative cancellation is that the anomalous induced 
by the QCD dynamics independent section of the first term of ()2.45|) was cancelled by the 
second term of ()2.45|) which is also QCD dynamics independent. As a consequence, one 
obtains the anomalous action from the QCD dynamics dependent section of the first term 
of f)2.45|) . In such a cancellation, the dependence of the Wess-Zumino term on the QCD 
can be seen explicitly. By considering the totally calculation given in the first subsection of 
section IIII[ we believe that the second cancellation is more reasonable. 

The QCD dynamics dependence of the Wess-Zumino term reflects its scale dependence. 
At low energy region, the dynamical mass is dominate and chiral symmetry is broken dy- 
namically. While at high energy region, the chiral symmetry is explicitly broken due to the 
small quark mass, and the chiral symmetry is preserved when the quark mass is neglected, 
then, there will be no anomalous meson process. This situation is very similar to the chiral 
anomaly studied in ref . [35] based on the loop regularization 3(| , in which both the massless 
and massive QCD will be anomaly free when the sliding energy scale fi s is large enough. 

In our calculation, we dealt with the fermion determinant explicitly by using the expan- 
sion of the QCD dependent functions A(k 2 ) and B(k 2 ) around the external sources. Besides 
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this method, the anomalous action from the term depending on the QCD dynamics can also 
be evaluated with Fujikawa's path integral method 34J, the anomalous action based on the 
general quark propagator can also be calculated with this method. But in this method, the 
behavior of the quark self-energy and the wave function renormalization constant should be 
determined under the limit r — > oo with r as the regulator introduced in this method, this 
makes the effects of the dynamical functions obscure. 

Generally speaking, the quark propagator with external sources is much more complex 
than that we considered here. Explicitly, it should be Lorentz and gauge covariant|31| and 
can be decomposed in the form of the Lorentz structure, i.e., the summation of the vector, 
axial- vector, scalar, pseudoscalar and tensor terms. It is noticed that the conclusion is much 
more complex and the physical meanings of the coefficients still deserve to be investigated, 
we shall discuss it elsewhere. 

Besides the anomalous action with pseudoscalar mesons, this method can be extended 
to the processes with resonances such as the vector mesons and axia-vector mesons are 
incorporated into the chiral effective theory by the hidden local symmetry method 21, S|- 
The leading gauge invariant anomalous action with external gauge fields was constructed 
infl by adopting the topological method with considering the t'Hooft matching condition. 
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